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P_i . We present optical {g' , Rc, and Ic) to near- infrared (J) simultaneous photometric 

observations for a primary transit of GJ3470b, the second lowest-mass planet among 
[ transiting planets around nearby M dwarfs, by using the 50-cm MITSuME telescope 

^P^' and the 188-cm telescope both at Okayama Astrophysical Observatory. From these 

2 ■ data, we derive the planetary mass, radius, and density as 14.1 ± 1.3 Mq, 4.321q^q 

c/3 ■ R^, and 0.94 it 0.12 g cm~'^, respectively, thus confirming the low density that was 

^ reported by Demory et al. based on the Spitzer /IRAC 4.5-//m photometry (0.72l;Q j^2 

g cm~^). Although the planetary radius is about 10% smaller than that reported 
^ i by Demory et al., this difference does not alter their conclusion that the planet pos- 

^ I sesses a hydrogen-rich envelope whose mass is approximately 10% of the planetary total 

I mass. On the other hand, we find that the planet-to-star radius ratio {Rp/Rg) in the 

!>: ■ J band (0.07577l[j-[}|][|^^) is smaller than that in the h (0.0802 ± 0.0013) and 4.5-^m 

O 
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(0.07806lJ];JjJ]J]^|) bands by 5.9±2.0% and 3.0± 1.2%, respectively. A plausible explana- 
tion for the differences is that the planetary atmospheric opacity varies with wavelength 
due to absorptions and/or scattering by atmospheric molecules. Although the signifi- 
cance of the observed Rp/Rs variations is still low, if confirmed, this fact would suggest 
that GJ3470b does not have a thick cloud layer in the atmosphere. This property would 
offer a wealth of opportunity for future transmission spectroscopic observations of this 
planet to search for certain molecular features without being prevented by clouds. 

Subject headings: planetary systems — planets and satellites: atmosphere — planets 
and satellites: individual(GJ3470b) — stars: individual(GJ3470) — techniques: photo- 
metric 



INTRODUCTION 



Transiting extrasolar planets provide much valuable information, not only planetary mass and 
radius, but also planetary atmospheric composition. Because the optical thickness of a planetary 
atmosphere varies with wavelength depending on the atmospheric composition, o ne can probe atmo- 
spheric constitueri ts by measuring transit radii with different wavelengths (e.g. ISeager &: Sasselov 



2000 



Brown 



2001 



So far, this technique, known as the transmission spectroscopy, has been applied for sev- 
eral transiting hot Jupiters orbiting nearby bright host stars, as represented by HD 209458b and 
HD 189733b. In the atmosph ere of HD 209458b, past observations detecte d m any absorption and 
scattering fe atures such as Na (ICharbonneau et al.ll2002l:ISnellen et al. I I2OO8I1. H dvidal-Madiar et al. 



20031). O, C (IVidal-Madiar et al 



2008b!), H2O (jBeaulieu et al.ll20ld ). and CO (jSnellen et al.ll20ld ). On the contrary, HD 189733b 



20041 ). TiO, V O (iDesert et al.ll2008l 1. H2 (jLecavelier Des Etanes et al 



regions (jPont et al. 



2008 



Sine et al. 


2009; 


Desert et al. 


2011b: 


Gibson et al. 


2012), while absorr 


were detected ( 


iledfield et al. 


20081; 


Lecavelier Des Etangs et al. 


2OI0I 



This featureless spectrum has been inter preted as Rayleigh scattering due to high-altitude haze 
which dominates over molecular features (jLecavelier Des Etangs et al.ll2008al ). 



Recently, it has become possible to expand this technique to low-mass planets, often referred 
to as exo-Neptunes (10 < Mp/M^ < 30) and super-Earths (Mp/M^ < 10), thanks to the discov- 
eries of transiting low-mass planets around nearby low-mass stars (M dwarfs). Because M dwarfs 
have smaller radii compared to Sun- like stars, they show deeper transits for a same-sized transiting 
planet, enabling us to obtain high signal-to- noise-ratio transit signals, which are need ed for trans- 



mission spectroscopy, ev en for small pl anets. The first two such examples are GJ436b (iButler et al 



20041 : iGillon et al.ll2007l ) and GJ1214b (jCharbonneau et al.ll2009l ). both of which are orbiting nearby 
M dwarfs. They can be thought to be the representatives of exo-Neptunes and super-Earths, re- 
spectively, in the sense of their masses (~23 and ~6.6 M^, respectively). As for GJ436b, no firm 
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molecular feature has been detect ed by transmission spectroscopy probably because of difficultie s 



due to instrumental systematics (jGibson et alj l201lh and stellar activity (IKnutson et al 



2011 



while a methane-poor and CO-rich atmospheric model has been supported (jKnutson et al.ll201ll ) 
by cor nbining the interpreta tion of the dayside spectrum obtained via secondary eclipse observa- 



tions ([Stevenson et al.ll20ld ). The super-Earth GJ1214b has received much more attention since 



its discovery. Intensive observations by many observational groups have revea l ed that GJ1214b 
has a fiat spectru i n over optical to infr a red wavelengths (e . 



Desert et al 



2nila 



Bean et al 



2011 



Berta et al.ll2012l : Ide Mooij et al.ll2012l : iNarita et al.ll2012l : iFraine et alj|2013l ). This flat spectrum 
has been interpreted as the consequence of either a water-dominated atmosphere or a hydrogen- 
dom inated but hazy /cloudy atmosphere, although which model is correct is still an open question 
(e.g. iHowe fc Burrowsll2012l ). 



GJ3470b, the target of this paper, i s the third low-mas s (< SOMe) transiting planet discov- 
ered among nearby (< 35 pc) M dwarfs (jBonfils et al.ll2012l ) that provides a great opportunity of 
extending the atmospheric study for low-mass planets. Because its mass, ~ 14M0, is intermediate 
between those of GJ436b and GJ1214b, this pl anet should be a good sample for comparative study 
of atmospheric properties of low-mass planets. iDemorv et al.l (j2013l ) (hereafter D13) have recently 
reported follow-up observations for this system including Spitzer /\R,KC 4.5-^m photometry of 
two primary transits of GJ3470b. They precisely determined the planetary mass and radius as 



13.9I^4M^ and 4.83_q2i-R®, respectively, reveahng its very low density {pp = 0.721q;^2 S cm~'^). 
This implies that GJ3470b has a light-gas atmosphere with enlarged spectrum features, thus mak- 
ing this planet as an attractive target for transmission spectroscopic studies. Nonetheless, no such 
study for this planet has been reported yet. 

Here we present optical-to-near-infrared simultaneous photometric observations for a primary 
transit of GJ3470b obtained by using the 188-cm telescope and the 50-cm MITSuME telescope both 
at Okayama Astrophysical Observatory. The simultaneous observations have a great advantage 
for transmission-spectroscopic study. If the host star has cool star-spots on its surface, transit 
depth (square of planet-star radius ratio) can vary with time due to the variations of apparent 
luminous area of the star according to the stellar rotation and appearing/ vanishing of the star- 
spots (e.g. iPont et al.ll2008l ). This changes the apparent planetary radius with time. Therefore, the 
simultaneous observations enable us to investigate the wavelength dependence of planetary radius 
without concern for this star-spot effect. 

The rest of this paper is organized as follows. We describe our observations in Section [21 
followed by data reduction and analysis in Section [3] and Section HI respectively. We discuss 
implications from our results in Section (Sj and summarize this study in Section [H 
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2. OBSERVATIONS 

2.1. ISLE J-band Observations 

We conducted J-band photometric observations for GJ3470 on the expected transit night 
of 2012 November 15, by us ing the near-infrared imaging and spectroscopic instrument ISLE 



(jYanagisawa et al.ll2006l . l2008l ). which is mounted on the Cassegrain focus of the 188 cm telescope 
at Okayama Astrophysical Observatory in Japan. ISLE has a 1,024 x 1,024 HgCdTe HAWAII-l 
array with the pixel scale of 0".245 pixel"^, providing the field of view of 4'. 3 on a side. For rel- 
ative photometry, we introduced a comparison star (2MASS:07591316+1525479, J=8.73) onto the 
detecter simultaneously with the target star GJ3470, where the two stars are separated by 2'. 9. 
At the beginning of the observations, the stellar positions on the detector were carefully set so 
that the both stellar images did not cover any bad pixels. During the observations, the telescope 
was defocused so that the full width at half maximum of the stellar point spread function (PSF) 
was ~20 pixel, or ~5", in order to mitigate the incompleteness of pixel-to-pixel-sensitivity (flat- 
fielding) correction and to extend the exposure time as long as possible while avoiding saturation. 
The exposure time was set to 30 sec, with which the peak analog-to-digital-unit (ADU) count of 
the brighter star (the comparison star) is well below the threshold of 25,000 ADU, above which 
the AD converting relation deviates more than 1% from a linear relation. The dead time including 
readout time for each exposure was 6 s. An observing log is shown in Table [TJ 

The 188-cm telescope is equipped with an offset guider system on the Cassegrain focus, which 
can correct the telescope tracking error every one second at minimum by using a fast-readout 
CCD camera which can keep track of a bright star coming into the surrounding area of the main 
detector. However, it had been recognized from past observations with ISLE that the stellar centroid 
positions on the ISLE detector slightly (a few pixels) changed over several hours even when the 
guiding system was activatec0. This slight stellar displacement had introduced systematic errors 
into photometry. In order to correct this effect, we have developed an offset-correcting system which 
corrects the reference point of the guide star on the guiding camera, by calculating the displacement 
of stellar centroids on the last ISLE images relative to a reference image. Each correction needs an 
additional dead time for ISLE of 8 s. For the transit observations of GJ3470b, we activated this 
system with the frequency of the offset correction as once in 10 ISLE exposures. As a result, the 
stellar centroid change during the observations was well suppressed with 1.1 and 0.6 pixels in root 
mean square (RMS) for the X (right ascension) and Y (decrination) directions, respectively. The 
centroid variation with time on the night is shown in Figure [H 



We attributes the cause of this stellar displacement to a relative mechanical offset between the ISLE detector 
and the guiding CCD camera depending on the telescope position, rather than other effects such as the differential 
atmospheric refraction between infrared and optical wavelengths. 
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2.2. MITSuME Optical Observations 

At the same time with the ISLE J-band observations, we also conducted photometric observa- 
tions for the same transit of G J3470b in optical wavelengths by using the 50-cm MITSuME telescope 
at Okayama Astrophysical Observatory. The telescope is equipped with three 1024 x 1024 pixels 



CCD cameras, enabling to obtain /c, Rc, and g' band images simultaneously (IKotani et al.ll2005l : 



Yanagisawa et al.ll2010l ). The each camera has the pixel scale of l".5 pixel ^ for the field of view 



of 26'x 26'. 

The telescope was defocused so that the FWHM of stellar PSF is about 10 pixels, or ~15", 
for the same purpose as the ISLE observations. We note that the contamination light from objects 
surrounding the target star is negligible because there is no object brighter than 20 and 16 magni- 
tudes in any band within 10'' and 2 0'', respectively, from GJ3470 in the SDSS photometric catalog 



(jAdelman-McCarthv &: et al.ll201ll ). The exposure time was set to 60 s. The dead time including 



readout time for each exposure was 3 s for all bands. The observing log is compiled in Tabled! 

Because this telescope has no mechanical auto-guiding system and previously caused large 
tracking error (~100 pixels) over several hours, we have developed a self-guiding software which 
calculates the displacement of stellar centroid positions on the last-observed /c-band image relative 
to a reference image, then feeds back it to the telescope to correct the tracking error soon after the 
last image is obtained. By using this self-guiding software, the stellar centroid displacement during 
the observations for GJ3470 were kept with 1.4 and 0.5 pixels in RMS for X (right ascension) and 
Y (declination) directions, respectively, for all the three bands. 



3. DATA REDUCTION 

3.1. Reduction and Baseline Correction for the ISLE data 

The obtained ISLE images are reduced by the standard procedure, including dark-image sub- 
traction and flat-fielding correction. For the flat-fielding correction, 100 dome-flat images obtained 
on the observing night are used to create a flat-fielding image. Then, aperture photo metry is done 



for th e target and comparison star on the reduced images by using a customized tool (jFukui et al 



20111 ). applying a constant aperture radius for all images. A light curve is created by dividing 
GJ3470's fluxes by comparison's ones. We eliminate apparent outliers due to such as passing 
clouds and cosmic-ray hitting from the light curve. The time for each data point is assigned as the 
mid time of each exposure in the time system of Barycentric Julian Day (BJD) based on Barycen- 
tric Dynamical Time (TDB), which is converted from the time stamp recorded on the FITS header 
in the t ime system of Julian Day (JD) based on Coordinated Universal Time (UTC) by using the 



code of lEastman et al.l (|20ld ) 



In order to select an appropriate aperture radius, we create a number of trial light curves with 
changing the aperture radius by 1 pixel, and evaluate the dispersion of the out-of-transit (OOT) 
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part of these trial light curves. We find that the light curve produced with the aperture radius of 
34 pixels gives the minimum dispersion, and therefore use this light curve for further analyses. The 
selected light curve is shown in Figure [TJ Subsequently, in order to correct the systematic trend in 
the light curve, which is apparent in the OOT light curve, we take the following procedure. First, 
we fit the OOT light curve with several baseline models given by the following formulae: 

Fbase = A;o X lO-O-^^™— , (1) 

Amcorr = '^hXi, (2) 
i=l 

where -Fbase is the baseline flux, {X} are variables, and {k} are coefficients. Generally, systematic 
trends could arise due to airmass change, slow variability of GJ3470 and/or the comparison star 
themselves, and stellar displacement on the detector. Therefore, for the variables {X}, we test 
several combinations of z, t, t^, dx, and dy, where z is airmass, t is time, dx and dy are the relative 
centroid posit ions in x and y directions, respectively. Next, we evaluate the Bayesian information 



criteria (BIG: ISchwarzjll978l ) for the respective baseline models. The BIG value for the OOT light 



curve is given by BIGoot = Xoot + iVoot , where Xoot the value of the baseline fit for the OOT 
part, k is the number of free parameters, and iVoot = 155 is the number of OOT data points. We find 
that the {X} = {z, dx, dy} model gives the minimum BIGoot value among all the baseline models, 
and therefore apply this model to correct the systematic trend in the entire light curve. In Table EJ 
we list the BIGoot values and the RMS of the OOT light curve (RMSoot) for the four representative 
baseline models of {X} = {z}, {z,dx,dy}, {z,t,dx,dy}, and {t,t'^,dx,dy}. Noticeably, including 
dx and dy in {X} provides a significant improvement on both the BIGoot and RMSoot values, from 
BIGoot = 224.2 and RMSoot = 1-61 x lO'^ for the {z} model to BIGoot = 171.4 and RMSoot 
= 1.35 X 10^'^ for the {z,dx,dy} model. This fact indicates that the stellar displacements on the 
detector produces significant systematics on the photometry, and the displacements during our 
observations were enough suppressed such that the systematics can be corrected, owing to the 
development of the offset-correction system (Section 12. Ih . In Figured! we also plot the the best-fit 
OOT light curve with the baseline model of {X} = {z, dx, dy}, along with the time variations of z, 
dx, and dy. The corrected light curve is shown in Figure [21 

The selected baseline model of {X} = {z, dx, dy} is, however, not so significantly favored 
compared to other baseline models of {t, z,dx,dy} and {t,t'^,dx,dy}, given the BIGoot difference 
of only ~ 5. Therefore, the choice of baseline model could be a cause of a systematic error on the 
final result, i.e., transit parameters. In order to evaluate the impact of this systematic effect, we 
fit the entire light curves corrected by different baseline models with a transit model, and compare 
the resultant parameters. The transit parameters we use are the planet-star radius ratio, Rp/Rg, 
the mid transit time, T^, the semi-major axis normalized by the stellar radius, a/Rg, the impact 
parameter, b = acosi/Rs, where i is the orbital inclination, and the coefficients for stellar limb- 
darkening effect. We adopt the orbital eccentricity and orbital period as e = and P = 3.33671 
days, respectively, from D13. For a stellar limb-darkening model, we use the quadratic limb- 
darkening raw, I{fi) = 1 — ui(l — fi) — ti2(l — ^)^, where / is the intensity, /i is the cosine of the 
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angle between the line of sight and the hne from the stellar center to the position of the stellar 
surface, and ui and U2 are the coefficients. Because ui and U2 are heavily correlated and can 
not independently be constrained well, we let ui free while fix U2 at a theoretical value during 
the fitting process; we adopt U2 = 0.255, witch is the mean of the two values for {T^s, loggs} = 
{3600, 4.5} and {3600, 5.0}, where Teff (K) is t he stellar eff e ctive temperature and log^s (cgs) is 
the stellar surface gravity, in the table given by I Claret et al.l (l2012l ) (The Tefj and log Qs values for 



GJ3470 were derived by D13 as 3600 ±100 and 4.658ifc0.030, resp ectively). For calculating a transit 
model, we use the an alytic formula given bv lOhta et al.l ( 20091 ). which is equivalent with that of 



Mandel &: Agoll (|2002l ) when using the q uadratic limb-dar kening raw. The best-fit parameters are 



determined by the AMOEBA algorithm ([Press et al.lll992l ). In Figure [21 the best-fit transit model 
for the light curve corrected with the {X} = {z, dx, dy} model is plotted. The photometric errors 
are rescaled so that the reduced value for the best-fit transit model becomes unity. 

The derived best-fit values and uncertainties of Rp/Rg, a/Rg, and T^, as well as the values 
are listed in Table [2j While the value for the {X} = {z} model is relatively large as expected, 
those for the other three models of {z,dx,dy}, {z,t,dx,dy}, and {t,t'^,dx,dy} are very close to 
each other. However, the Rp/Rg values for these three models are slightly different, with the 
largest difference of 0.00081 between {z,dx,dy} and {t,t'^,dx,dy}, which is comparable to the 1-a 
uncertainties of Rp/Rg- Therefore, we consider this discrepancy as a systematic error in Rp j Rg , 
and will be taken into account in the study of wavelength dependency of Rp/Rg discussed in Section 
15. 2i On the other hand, the baseline-model dependences ona/Rg, b, and Tc are negligible compared 
to their 1-a uncertainties. 



3.2. Reduction and Baseline Correction for the MITSuME data 

The obtained MITSuME images are reduced in the same way as the ISLE images. For the 
flat-fielding correction, 27 twilight-flat images for each band obtained on the observing night are 
used to create each flat-fielding image. We then apply aperture photometry for GJ3470 and dozens 
of bright stars on the reduced images, with a number of trial aperture radii incremented by 0.5 
pixels. Among a number of combinations of aperture radii and comparison stars, we first visually 
select several good combinations which produce less-dispersed light curves with respect to a transit 
signature. Next, we fit each trial light curve including the transit part with a transit + a tentative- 
baseline model, and select the least-dispersed light curve as the most appropriate one for each band. 
The reason of including the transit part is that, unlike the J-band data, the MITSuME data do 
not have enough number of OOT data. For the tentative baseline model, we use Equation (1) and 
(2) with {X} = {z}. As for the transit model, we use the same parameterization as Section [3.11 
For the Iq- and Rc-hand data, we fix U2 at the theoretical values of 0.338 and 0.322, respectively, 
while let ui free during the fitting process. For the ^('-band data, we fix both ui and U2 at the 
theoretical values of 0.486 and 0.289, respectively, because it is too poor to constrain ui nor U2 by 
the current data. These theoretical values are derived by the same way as for the J-band data. In 
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addition, we force a/Rs and h to the values derived from D13, namely 13.4 and 0.40, respectively. 
This is because the Spitzer^s data used by D13 is more precise than the MITSuME's data, and we 
are more focusing on studying the wavelength dependence of Rp/Rg (see Section [5.2p rather than 
constraining other parameters. As a result, we select the light curves produced with the aperture 
radii of 7.0, 8.0, and 6.0 pixels, and the number of comparison stars of 4, 10, and 8, for /(., Re, and 
g' bands, respectively. We note that these results are robust over different baseline models. 

Subsequently, in order to select appropriate baseline models for correcting systematic trends, 
we again fit the selected light curves with a transit + several baseline models. At this time, we test 
the three baseline models of {X} = {z}, {z,t}, and We do not include dx and dy in {X} 

because we do not find any improvement for any bands when we preliminary fit the light curves by 
including dx and dy. In Table O we list the resultant BIC and RMS values for the three baseline 
models for each band. For the Iq and g' bands, we find that the baseline model of {X} = {z} gives 
the minimum BIC values, and therefore adopt it for correcting the systematic trends in these light 
curves. On the other hand, for the iic-band data, all the three baseline models give almost the 
same BIC values, indicating that the three models have the same statistical significance. So, we 
decide to select the model {X} = {z, t} among the three models, based on the fact that the RMS 
values for {z,t} and {t,t^} (2.46 x 10^'^) are slightly better than that for {z} (2.50 x 10~^), and it 
is physically more straightforward to include airmass in the variables rather than just a polynomial 
function of time. 

As also discussed in Section ISTTl the selection of baseline models could affect on the final result. 
In Table El we also list the resultant Rp/Rg values for the respective models for each band. The 
largest differences oi Rp/Rg are 0.0014, 0.0060, and 0.009, for the Ici Rc, and g' bands, respectively. 
The differences for the Ic and g' bands are comparable to and less than the l-cx uncertainties of 
Rp/Rg, respectively, while that for the Rc band is more than 3 times the 1-a uncertainty of Rp/Rg. 
These possible systematic offsets will be considered in the discussion on the wavelength dependence 
of Rp/ Rg in Section 15.21 

In Figure El we show the uncorrected light curves and the best-fit transit+baseline models of 
the three bands, as well as their residuals. After the baseline for each light curve is corrected, error 
bars are rescaled so that the reduced ^oi a transit-model fit becomes unity. 



ANALYSIS 



In order to derive the final transit parameters and their uncertainties, we analyze the corrected 
J-, Ic-, Rq-, and 5f'-band light curves simultaneously by the Markov Chain Monte Carlo (MCMC) 



method by using a customized code (jNarita et al.ll2012l ). In this analysis, we treat b, a/ Rg, and Tc 
as common parameters for all the four light curves, while Rp/Rg, ui, U2 as independent parameters 
for the respective light curves so as to take the wavelength dependences of these parameters into 
account. The values of b, a/Rg, Tc, and each Rp/Rg for each band are left free during the MCMC 
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process. The ui and U2 values are treated the same way as previously; ui for the J, /c, and Rc 
bands is left free, while ui for the g' band and U2 for all the bands are fixed at the theoretical 
values. We adopt e = and P = 3.33671 days in the same way as previously. 

In a MCMC chain, we start with a set of parameters which provides the minimum value 
determined by the AMOEBA algorithm. The value is given by 

2 f fohsA /model \ /„n 

^=^1 a. J' 

where /obs,i is the i-th photometric flux, is the i-th photometric error, and /model is the model 
flux calculated from the given parameter set. Then, the next parameter set is randomly selected, 
where we assume gaussian distributions whose means are the respective current values and standard 
deviations are proportional to their 1-a errors. This parameter set is always accepted when it gives 
Ax^ < 0, where Ax^ is the difference between the next parameter set and current one; otherwise 
it is accepted according to an acceptance probability which is given by p = exp^-^^/^. This process 
(link) is repeated for 10^ times in a chain, from which the first 10^ links are removed as a burn- 
in portion. The widths of the gaussian distributions for jumping parameters are set so that the 
acceptance ratio becomes about 25%. We run 10 independent chains, and create merged posterior 
distributions of the respective parameters. We define l-a statistical uncertainties as the range of 
parameters between 15.87% and 84.13% of the posterior distributions. The resultant values and 
uncertainties are listed in Table IH and two-dimensional 68.3% and 99.7% confidence regions for 
selected parameters drawn from the posterior distributions are shown as red contours in Figure 
m The derived model light curves for the respective bands are shown in Figure [U along with the 
observed data binned in 5-minute intervals. 

We confirm that the values of b, a/Rs, and Rp/Rs for Jc, Rc, and g' bands are all consistent 
with those reported by D13 within 2 a. On the other hand, we find that the Rp/Rg value for the 
J-band data, 0.07536 ± 0.00079, is inconsistent with that for the Ic-band data of 0.0797 ± 0.0014 
by 2.6 a, and also inconsistent with that for the Spitzer^s 4.5-^m data of 0.07806jiQ'QQQ54 by 2.8 
a. One possibility of these discrepancies is a systematic effect originated from stellar activity (e.g. 



Czesla et al.l l2009l ). Occulting cool star-spots by a transiting planet can induce a bump in the 
transit light curve, easily affecting the fitted parameters. Although we cannot identify any such 
feature in any of our light curves nor the Spitzer^s ones, if the planet occulted cool spots existing 
near the stellar limb, they would be difficult to be identified, and cause systematic offsets on the 
derived transit parameters such as b, which is weakly correlated with Rp/Rg (see Figured]). This 
might be the case, because the b value derived from our data (< 0.281 as 1-a upper limit) is slightly 
inconsistent with that by D13 (0.40lo;og) by more than 1 o". In order to reduce such a systematical 
effect on Rp/Rs, we repeat the MCMC analysis by using the b value reported by D 13 as a prior 
information. In this analysis, we use the following function instead of Equation ([3|) for calculating 
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the value: 

2 I f fohsA ~ /model \ , f ^ ~ &prior \ | / . x 

where we adopt ftprior = 0.40, and (yb^^iar ~ ^-^^ ~ ^prior) > 0; otherwise crfep^..^^ = 0.08. The 
resultant values and uncertainties are summarized in Table HJ and the two-dimensional correlation 
maps for selected parameters are plotted as blue contours in Figure [H As a result, we find that our 
b and a/Rg values become consistent with those by D13 within 1 a. However, although the Rp/Rg 
values for our data (J, Ic, Rc, and g' bands) become to be closer to that for the 4.5-/im band, there 
are still disagreement between J and Ic by 2.5 a, and between J and 4.5 ^m by 2.6 a. In principle, 
these discrepancy could be due to the wavelength dependence of the planetary atmospheric opacity. 
This will be discussed in Section 15.21 



5. DISCUSSIONS 

5.1. Physical Parameters of the Planetary System 

In this section, we focus on the physical parameters of the planetary system. One of our goals 
in this section is to test the very low density of GJ3470b suggested by D13 with our independent 
observations, and therefore we here discuss along the results of the MCMC analysis without b prior. 

First, we refine the orbital period by a linearly fit to our Tc value and the DlS's one (they 
reported one reference transit mid-time). The resultant orbital period is P = 3.336651 it 0.000005 
days, which is consistent with P = 3.33671 ±0.00005 days deriv ed by D13 based on the two transits 



observed by Spitzer and four transits from (jBonfils et al.ll2012l ). with the difference of only 5. lit 4.3 



sec. This means that no transit timing variation exceeding 4.0 it 3.4 min is observed between the 
two epochs, which are separated by 47 transit epochs. 

Next, we investigate the stellar mass and radius. The stellar density, ps, can directly be 
derived from a/Rs and P via the following relation assuming a circular orbit: ps = 0.01342 x 
(a/i?s)^/(-P(days))^ Pq. From our MCMC results, we derive the stellar density as ps = 3.83!1^q32 
Pq, which is broadly consistent with the value derived by D13 (2.91^033 Pq) within 2 a, although 
our value is 32% larger than that by D13. The stellar mass and radius can be derived from the stellar 
density, combined with one or more additional information such as a stellar color (or temperature), 
luminosity (i.e., apparent brightness and distance to the star), and/or a stellar mass-radius relation. 
Because a trigonometric parallax for GJ3470 has not been measured so far, D13 solved the stellar 
mass, radius, and parallax si multaneously from p ,^ and VJHKg magnitudes, via two empirical 



relat ions of mass-luminosity (jPelfosse et al.l I2OOOI ) and diameter-color-luminosity (jKervella et al 



20041 ) relations. As a result, they derived the stellar mass and radius as Ms = 0.539lo;o43 M© and 



Rs = 0.568^gQ3[ Rq, respectively. However, the latter relation they adopted was derived from dwarf 
samples which contained relatively small number (< 15) of low-mass (K-M) dwarfs. In addition. 
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metallicity dep e ndenc e of this relation was not found in the samples (jKervella et al.ll2004l ). however, 
Bovaiian et al.l (|2012l ) have recently found, based on interferometric radius measurements for 33 
single K-M dwarfs, that a color-radius relation for low-mass stars clearly depends on metallicity. 
Therefore, the relation adopted by D13 could potentially h a,ve a systematic offse t. In Figure O 
we plot a mass-radius relation for low-mass stars provided bv lBovaiian et al.l (|2012l ) (green shaded 
region), which was derived based o n the radius measur e ment s for the 33 K-M dwarfs and an 
empirical mas s -lumi nosity relation of iHenry &: McCarthy! (119931 ) . which is consistent with that of 
Delfosse et al.l ([20001) • we also plot a mass-radius relation for GJ3470 drawn by connecting the two 
empirical relations that D13 adopted (light-blue shaded region), where Kg = 7.989 is used for the 
mass-luminosity calibration and H = 8.206 and {V — H) = 4.124 are used for the diameter-color- 
luminosity calibration. These two relations do not cross each other, ir nplying that sorne sys tematic 
offset could exist, possibly in the diameter-color-luminosity relation of lKervella et al.l (j2004l ) due to 
the small number of low-mass star samples and the unconsidered metallicity dependence. 



Therefore, instead of using the same 



way 



as D13, we solve Mg and Rg from ps by using 
the mass-radius relation of lBoyajian et al.l (120121 ) (E quation (4) in the pap er). We note that any 
metallicity dependence of this relation was detected (|Bovaiian et al.ll2012l ). In order to properly 
estimate M^, Rg and their uncertainties, we create probability distributions of these parameters 
by using a Monte Carlo (MC) technique. That is, the probability distributions are created by 
repeating the process that the values of a/Rg and three coefficients of the mass-radius relation are 
taken according to their respective probability distributions, and then Mg and Rg are solved. For 
the probability distribution a/Rg, the posterior distribution created from the MCMC analysis is 
used, while for those of the three coefficients, gaussian distributions whose standard deviations are 
set to their 1-a errors are used. The resultant median values and the 68.3% confidence intervals are 
calculated as Mg = 0-557toml Mq and Rg = 0.526to!^ll Rq. In Figure El we also plot the derived 
Mg and Rg values (red circle) as well as those reported by D13 (brown triangle). Our Mg and Rg 
values are both consistent with those by D13 within 1 a. However, this agreement is by coincidence, 
because the derived stellar density is slightly inconsistent between this work and D13 by more than 
1 a, and the adopted calibration methods for deriving stellar mass and radius are different and 
also inconsistent between the two. Thus, further follow-up observations to confirm the stellar mass 
and radius would be valuable; photometric transit observations could test the stellar density, and 
astrometric observations for measuring the trigonometric parallax would provide a new insight into 
the stellar mass/radius calibration. 

Finally, using the estimated Mg and Rg, we derive the relevant stellar and planetary parameters. 
For the planetary radius, we adopt the Rp/Rg value for J band. This is because not only the J- 
band light curve has the highest photometric precision among the four light curves, but also, as 
will be discussed in Section 15.21 the planetary radii for the optical bands could be enlarged due 
to putative atmospheric haze, and the smaller J-band planetary radius could represent a haze-free 
planetary radius. For the planetary mass, we use the following relation assuming a circular orbit: 
K' = {2'nGf/^Mps\ni/{Mg + Mpf'^, where K' is a parameter derived from radial-velocity data 
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and K' = 13.4 it 1.2 m d^/^ is adopted from D13. The Mp and Rp values as well as other stellar 
and planetary parameters (the stellar surface gravity, loggs, the semi-major axis, a, the planetary 
density, pp, and the planetary surface gravity, gp) are derived by the same MC fashion as previously. 
The derived values and uncertainties are listed in Table [5j The planetary mass, radius, and density 
are derived as Mp = 14.1 ± 1.3 Mj^, Rp = 4.32l[J;^o i^gj, and pp = 0.94 ± 0.12 g cm"!, respectively. 
Although the planetary density is ~30% larger than that by D13 (0.72]']q|2 g cm~^), it is still well 
lower than that of Uranus (1.27 g cm~^) de spite of its si r nilar mass (14.54 M^), confirming the 
low density of GJ3470b, firstly suggested by iBonfils et al.l (j2012l ) and recently uncovered by D13. 
We note that the ~30% difference of the planetary density comes mainly from the ~10% smaller 
planetary radius (the ~7% smaller stellar radius and the ~3% smaller planet-star radius ratio) 
compared to D13. In Tabled we also list the results derived from the MCMC analysis with h prior, 
for reference. 

In order to investigate the impact of the 10% difference of the planetary radius on the plane- 
tary bulk composition, we have modeled the internal structure of GJ3470b that is consistent with 
the planetary mass and radius derived in this study. To do so, we have assumed that the planet is 
composed of three layers, namely, a cloud/haze- free hydrogen-rich atmosphere with solar composi- 
tion, on top of an ice/water mantle, on top of a rocky core (the water/rock mass ratio is set to 3). 
Two-layered structure (hydrogen-rich atmosphere -|- r ocky core) is also ex amined. The equations 



of sta te adopted are SCvH EOS for hydrogen/heli um (jSaumon et al.lll995l). R-EQS (IFrench et al 



20091 ^ and SESAME EOS (ILvon Johnsonlll992l ') for water, and IValencia et al.l ((20071) for rock 
(details are described in Kurosaki et al. in prep). We have found that the mass fraction of the 
hydrogen-rich atmosphere ranges 5% to 9% in the three-layer models and 12% to 19% in the two- 
layer models. These values are similar with those derived in D13. Therefore, the 10% difference 
in radius does not alter their conclusion that the planet possesses a hydrogen-rich envelope the 
mass of which accounts for approximately 10% of the planetary total mass. We note that this 
conclusion may be consistent with recent theoretical prediction: If the atmosphere embedded in a 
protoplanetary disk grows in r nass beyond 10% of the planetary total mass, its accretion tends to 
proceed in a runaway fashion (jlkoma &: Horill2012l ). Thus, nebular-origin atmospheres of >10% of 
planetary mass are rarely detected. 



5.2. Wavelength Dependence of the Planet-star Radius Ratio 

In this section, we discuss the wavelength dependence of Rp/Rs based on the results of the 
MCMC analysis with b prior. In Figure El we plot the Rp/Rs values in the g' , Rc, Iq, and J bands, 
and that in the Spitzer^s 4.5-/im band reported by D13, along with the transmission curves for the 
respective passbands. The Rp/Rs values in the Ic and 4.5-/xm bands are larger than that in the J 
band by 5.9 it 2.0% and 3.0 it 1.2%, respectively. Those in the g' and Rc bands are also larger than 
that in the J band, although the uncertainties in these bands are relatively large. 

As discussed in Section [3 . 1 1 and [3 . 2 1 the choice of the baseline model could cause a systematic 
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offset in Rp/Rg. However, the alternative choices of the basehne models for the i?c-j Iq--, and J-band 
light curves work as expanding the discrepancies, as indicated by gray circles in Figure [HI meaning 
that the systematics due to the baseline selections cannot explain the observed discrepancies of 
Rp/Rs- 

Unocculted star-spots can also vary apparent Rp/Rs with time, due to the variations of ap- 
parent luminou s area of the star according to the stellar rotation and appearing/ vanishing of the 



star-spots (e.g. iPont et al.l 120081 ). Since the observational epochs are different with ~5 months 
apart, there is a possibility that the observed Rp/Rs difference between the J and 4.5-/im bands 
is caused by the star-spot effect. When a spotted host star induces the peak-to-valley variability 
of r = Rp/Rs, Ar, due to the unocculted star-spot effect, the host star itself will show a peak-to- 
valley flux variability of (1 + (Ar/r))^ — 1 with a period consistent with the stellar rotation period. 
Therefore, if the 3.0 ± 1.2 % difference of Rp/Rs between the J and 4.5-//m bands were caused by 
star-spots, the host star would have to show at least 6.0 it 2.3% peak-to-valley intrinsic variabil- 
ity. In order to estimate the maximum intrinsic variability of the host star, we have conducted a 
long-term photometric monitoring for GJ3470 by using the 50-cm MITSuME telescope, spanning 
~60 days. The resultant Ic-band light curve is shown in Figure [71 along with the sky-flux varia- 
tions. We have found that the RMS of the nightly-averaged fluxes is 0.33% and the peak-to-valley 
variability is merely ~1%, indicating that GJ3470 is not a very active star and star-spot variations 
cannot account for the observed Rp/Rs difference between the J and 4.5-;um bands, unless the stel- 
lar variability is unusually large only in infrared wavelength. In addition, the difference of Rp/Rs 
between the Ic and J bands cannot intrinsically be explained by the spot-rotation effect, because 
these data were obtained simultaneously. We note that a possible ~25-day periodic variation seen 
in the long-monitoring light curve could potentially be due to a star-spot rotation according to the 
stellar rotation, however, we have found that a part of the flux variations (JD > 2,456,258) could 
be correlated with the sky- flux variations, and therefore further monitoring is needed to confirm it 
(more will be investigated in a future paper). 

More likely scenario would be that the planetary atmospheric opacity varies with wavelength 
due to absorptions and/or scattering by atmospheric molecules. The relative change in tran- 
sit depth due to molecular absorptions as a function of wavelength can be approximated by 
AD{X) ~ 2nu{X)HRp/ R^, where n^( X) is the scal e factor depending on the molecular opacity, 
which can be ~10 for strong absorbers ( BrownlboOll ). and H is the atmospheric scale height given 



hy H = kT/^Qp, where k is Boltzmann's constant, T is the atmospheric temperature, and // is 
the mean molecular weight. For GJ3470b, H~250-400 km when assuming T = 500-650K and a 
solar-composition atmosphere (/i = 2.36 atomic mass). Therefore, AD(X) can be up to ~0.16%, 
which corresponds to the relative change in Rp/Rs of ~13%, indicating that the observed differ- 
ences of Rp/Rs can reasonably be explained by molecular absorptions. Indeed, there are strong 
absorption bands of CO and CO2 around 4.5;um, which makes Rp/Rs in the 4.5-/im band be larger. 
Therefore, the observational result that Rp/Rs in the 4.5 fim band is larger than that in the J 
band is qualitatively consistent with a CO- and/or C02-containing atmosphere. On one hand, the 
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larger Rp/Rg values in the optical {Iq, Rc, and g') bands compared to that in the J band could be 
explained by Rayleigh scattering due to mol ecular hydrogen or small-sized (<0 -l /^m) haze in the 
atmosphere , as is the case for HD 209458b (jLecavelier Pes Etangs et al.ll2008bl ) and HD 189733b 
(|Pont et al.ll2008l ). Especially, the transition wavelength below which the spectr um begins to rise 
due to Rayleigh scattering is around l//m for ~0.1-;um sized haze particles (e.g. iHowe &: Burrows 
2012I ). and therefore the observed Rp/Rs difference between the Ic and J bands could be explained 
well by such a hazy atmospheric model. In Figure El we also plot a model spectrum which qual- 
itatively fits the data (b lue dotted line). For the model spectrum, we adopt a model provided by 
Howe &: Burrowa (|2012l ) in which a 0.3 x solar-abundance atmosphere with T = 700K containing 
0.1-^m tholin particles with a density of 1000 cm~^ is assumed, and scale it so as to fit the scale 
height of GJ3470b. We note that a quantitative discussion and fine-tuning of this model is beyond 
the scope of this paper, and to do so, further observational data are required. 

The observed Rp/Rs variations also indicate that the planet would have no thick clouds in 
the sky, because if the atmosphere of GJ3470b were covered by thick clouds at hi gh altitudes 



a mu ch flat spectrum over optical to infrared wavelengths would be observed (e.g. iBerta et al 



2OI2I ). If true, this fact would offer a wealth of opportunity to probe many molecular features in 
the GJ3470b's atmosphere through high-precision transmission spectroscopic observations, without 
being prevented by thick clouds. On the other hand, a constant fit to the observed five Rp/Rg 
values yields the value of 11.6, meaning that the statistical probability of exceeding this 
value assuming gaussian errors is 2.1%. Thus, a fully-clouded atmospheric model still cannot be 
ruled out by the current data. 

We should note that this is the first report on transit observations by using the 188-cm 
telescope/ISLE instrument and the 50-cm MITSuME telescope. Our observations for GJ3470b 
demonstrate that these telescopes/instruments are quite useful for studying transiting planets, es- 
pecially for probing planetary atmospheres by simultaneous observations through the optical to 
near-infrared wavelengths. Nevertheless, the capability of the 50-cm MITSuME telescope is lim- 
ited due to the limited aperture size. It would be a great capability if 2-4 m class telescopes 
like the 188-cm telescope are equipped with a transit-dedicated camera which can take multi- 
color images through optical to infrared simultaneously. Th e GROND instrument, mounted on 
the MPG/ESO 2.2-m telescope in Chile (iGreiner et al.l 120081 ). is a pioneer for such an ambitious 
instrument, which can obtain seven images through g' to K bands simultaneously for the southern 
hemisphere. R ecently, simultaneous transit observations by using the GROND instrument have 
been reported (jMancini et al.l 120131 : Ide Mooij et al.l |2012| ) , demonstrating the usefulness of such a 
multi-color imager for transit observations. However, this instrument had originally been developed 
for catc hing gamma-ra.y-burs t afterglows and not been designed for high-precision transit obser- 
vations ([Mancini et al.l 120131 ). Future developments of such instruments, but more dedicated to 
transit observations, for 2-4 m class telescopes in northern hemisphere will provide fruitful results 
on the exoplanetary atmospheric studies. 
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6. SUMMARY 



We have presented optical (g' , Rc, and Ic) to near-infrared (J) simultaneous photometric 
observations for a primary transit of the hot Neptune GJ3470b, by using the 188-cm telescope/ISLE 
instrument and the 50-cm MITSuME telescope both at Okayama Astrophysical O bservatory. We 



have found that the planetary density is broadly consistent with that reported by iDemory et al 



(|2013l ) (D13), who measured it based on the Spitzer /IRAC 4.5-//m photometry, confirming its low 
density. Although the derived planetary radius is about 10% smaller than that reported by D13, 
this difference does not alter their conclusion that the planet possesses a hydrogen-rich envelope 
whose mass is approximately 10% of the planetary total mass. On the other hand, we have found 
that the planet-to-star radius ratio (Rp/Rs) in the J band is smaller than that in the Ic and 4.5- 
^m bands by 5.9 ± 2.0% and 3.0 ± 1.2%, respectively. These discrepancies cannot be explained 
by systematic effects due to baseline corrections. In addition, we have found from a long-term 
flux monitoring for GJ3470b that the intrinsic peak-to- valley stellar variability in Ic band is mealy 
~1%, indicating that the unocculted star-spot effect is unlikely to account for the observed Rp/Rs 
difference between the J and 4.5-/i bands. Instead, Rayleigh scattering due to molecular hydrogen 
or small-sized (^0.1 /zm) haze in the atmosphere and molecular absorptions such as CO and CO2 
could reasonably explain the observed Rp/Rs variations. Although the significance is still low, if 
these Rp/Rs variations are confirmed, this fact would suggest that GJ3470b would not have a thick 
cloud layer in the atmosphere, offering a wealth of opportunity for future transmission-spectroscopic 
observations to probe many molecular features in the atmosphere of the hot Neptune. 

We note that this is the first report on transit observations by using the 188-cm telescope/ISLE 
instrument and the 50-cm MITSuME telescope. Our observations demonstrate that optical-to-near- 
infrared simultaneous observations using such as these instruments are very useful for planetary 
atmospheric studies. Future developments of multicolor imagers for 2-4 m class telescopes in north- 
ern hemisphere would provide much more fruitful results for this field. 
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Fig. 1. — (Upper panel) The J-band uncorrected light curve of GJ3470b. The best-fit baseline 
model, which is as a function of airmass and stellar centroid displacements in X and Y directions, 
is over plotted as solid line in the GOT part. (Second panel) The stellar centroid displacements in 
X direction. (Third panel) The same as the second panel but in Y direction. (Bottom panel) The 
airmass variations. 
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Fig. 2. — (Upper panel) The baseline-corrected J-band light curve. The best-fit transit model 
derived in Section [3T] is shown as solid line. (Bottom panel) The residual light curve. The RMS of 
the residuals is 1.18 x 10^^. 
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Fig. 3. — (Left panel) The uncorrected light curves of MITSuME /c, Rc, and g' bands, from top to 
bottom. The best-fit baseline models are over plotted as solid lines. For display, -0.02 and -0.048 
are added in the Rc and g' light curves, respectively. (Right panel) The same as the left panel but 
residual hght curves. The RMS of these residuals are 1.78 x 10~^, 2.45 x 10~^, and 5.27 x 10~^, 
for /c, Rc, and g' bands, respectively. 
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Fig. 4. — The two-dimensional 68.3% (thick hnes) and 99.7% (thin hnes) confidence regions for 
selected parameters calculated from the posterior probability distributions as a consequence of the 
MCMC analysis. The gray (red in online version) contours show the results from the MCMC 
analysis without any prior, while the black (blue in online version) contours indicate those with a 
b prior. 
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Fig. 5. — The overall light curves after the baseline correction, the error normalization, and 5-min 
binning. The J-, Iq-, Rc-, and i^'-band light curves are shown from top to bottom. The offsets of 
-0.08, -0.18, and -0.28 are added for the Ic-, Rc-, and ^('-band light curves, respectively, for clarify. 
The best-estimated transit models derived from the MCMC analysis without b prior are shown as 
solid lines. 
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Fig. 6. — Comparison of Mg and Rs calibrations between this work and D13. The stehar density 
derived in this work (MCMC analysis without b prior) and that by D13 are indicated by lower 
(red in online) and upper (brown in online) lines. Solid and dashed lines represent median and 1 a 



uncer tainties, respectively. The mass-radius relation we adopted (Equation (4) in iBoyajian et al 



(|2012l )) is shown as right (green in online) shaded regio n (1-a region), while the mass-radius relation 
for G J3470 drawn by connectin g a mass-luminosity (jPelfosse et al.l 12000 ) and a diameter-color- 
luminosity (iKervella et al.l l2004l ) relations are shown as left (light-blue in online) shaded region, 
where 2% uncertainty in Rs is adopted. The Mg and Rs values derived in this work and by D13 
are indicated as circle and triangle, respectively. 
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Fig. 7. — (Upper panel) A long-term Jc-band light curve for GJ3470 obtained by the 50-cm MIT- 
SuME telescope. Transit part is eliminated in this light curve. Individual 60-sec exposure data and 
nightly-averaged data are indicated as gray and black points, respectively. RMS of the black points 
is 0.33%, while peak-to-valley variations is The ±0.5% variation level is shown as dashed 

line. The transit-observed date is indicated by an arrow. (Bottom panel) Sky-flux variations. A 
possible correlation between the light curve and sky flux can be seen after JD-2450000 ~ 6258. 
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Fig. 8. — (Top panel) Wavelength dependence of the observed Rp/Rg values for GJ3470b. The data 
derived from the MCMC analysis with b prior are indicated as filled circles, with green, yellow, 
orange, and red for the g', Rc, Ic, and J bands, respectively. The Spitzer /IRAC 4.5-/xm data 
obtained by D13 is shown as brown triangle. The horizontal values and error bars denote the 
weighted center and the width at half maximum of the respective transmission curves shown in 
the bottom panel, respectively. The horizontal dashed line indicates the weighted mean of the five 
Rp/Rs values (=0.07761), representing an approximated fully-clouded atmospheric spectrum. The 
blue dotted line indicates a model spectrum for a 0.3 x Solar-abundance cloud-free atmosphere with 
T = 700 K and containing 0.1-^m sized t holin particles with a den sity of 1000 cm~^, being drawn 
by scaling a model spectrum provided by iHowe &: Burrows! (j2012l ) so as to fit the scale height of 
GJ3470b. The blue squares indicate the integrated cloud-free model spectrum over the respective 
passbands. The gray circles indicate the impact of possible systematics due to the baseline- model 
selection; the original value is shifted by the difference of Rp/Rs that would be produced if an 
alternative baseline model was adopted. (Bottom panel) The transmission curves for the respective 
bands. 
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Table 1. Observing Log 



Observing Date Filter Telescope Exp. time [s] A^aii ^ Airmass variation 



2012 November 15 J 

2012 November 15 h 

2012 November 15 Rc 

2012 November 15 g' 



OAO 188 cm 30 

MITSuME 50 cm 60 

MITSuME 50 cm 60 

MITSuME 50 cm 60 



352 1.42 1.06 ^ 1.12 

165 1.28 1.06 1.11 

162 1.28 1.06 1.11 

164 1.28 1.06 1.11 



The total number of observed data used for analyses, i.e., after omitting outliers. 



Table 2. Fitting results for J-band data for different baseline models ^ 



Variables 
{X} 


BICoot 


RMSoot 
(xlO-^) 


Xall 


Rp / Rs 
(xlO-2) 


a/Rs " 


b " 


[BJDtdb-2456247] 


z 


224.2 


1.61 


449.3 


7.754 ± 0.098 


14.91 tl% 


0.00 ± 0.21 


0.29949 ± 0.00025 


z, dx, dy 


171.4 


1.35 


299.5 


7.513 ± 0.082 




0.00 ± 0.24 


0.29971 ± 0.00022 


z, t, dx, dy 


176.3 


1.35 


299.5 


7.483 ± 0.082 


14.87 


0.00 ± 0.24 


0.29974 ± 0.00022 


t,t^,dx,dy 


176.4 


1.35 


299.9 


7.432 ± 0.083 


14.88 t°;ll 


0.00 ± 0.24 


0.29973 ± 0.00022 



The fitting results for tfic adopted baseline model are indicated as bold text. 

^ The subscript 'all' denotes that the entire light curve is used for the calculation. 

The uncertainties are calculated after the flux errors are rescaled so that the reduced Xaii value is unity. 

The b value is allowed to be minus in the fitting process, but the absolute value is used for modeling a transit 
light curve. 
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Table 3. Fitting results for the MITSuME 
data for different baseline models ^ 



Variables BIG RMS Rp/Rs 
{X} (xlO-3) (xlO-2) 



MITSuME 7c 



z 


185.5 


1.81 


7.97 ± 0.13 




192.1 


1.82 


8.06 ±0.12 


t,t^ 


191.8 


1.82 


8.11 ±0.12 




MITSuME Rc 




z 


186.5 


2.50 


8.32 ±0.17 


z,t 


186.5 


2.46 


7.72 ± 0.17 


t,t^ 


186.5 


2.46 


7.74 ±0.17 




MITSuME g' 




z 


184.5 


5.12 


8.16 ± 0.30 




189.6 


5.12 


8.08 ± 0.30 




189.7 


5.12 


8.07 ± 0.30 



The fitting results for the adopted baseline models 
are indicated as bold text. 
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Table 4. MCMC results for transit parameters 



Parameter 



D13 



This work (w/o b prior) This work (w/ b prior) 



Impact parameter, b {= a cos i/Rs) 
Scaled semi-major axis, a/Rs 
Mid transit time, Tc [BJDtdb - 2450000] 
Planet-to-star radius ratio 
Rp/Rs (4.5Atrn) 
Rp/Rs (J) 
Rp/Rs {Ic) 
Rp/Rs (Re) 
Rp/Rs (g) 

Limb-darkening coefficients 
Ml (J) 
U2 (J) 
Wl {Ic) 
U2 {Ic) 
Ul {Rc) 
U2 {Rc) 

Ul {g') 
U2 {g') 



0.40 



+0.06 



+0.55 



13.42 

6090.47690 ± 0.00015 



0.07806t°;rg 



< 0.281 (1-cr upper limit) 
14.70 

6247.29954 ± 0.00019 



0.07536 ± 0.00079 
0.0797 ± 0.0014 

0.0770 ± 0.0019 
0.0800 ± 0.0030 

0.149 ± 0.074 
0.255 (fixed) 
0.20 ± 0.11 

0.338 (fixed) 
0.26 ± 0.15 
0.322 (fixed) 
0.486 (fixed) 
0.289 (fixed) 



n QQ7+0.067 
U.OO(_g 070 



-0.39 

6247.29951 ± 0.00020 



n 07577 +0-00072 

U.U(0( ( -0.00075 

0.0802 ± 0.0013 

0.0776 ± 0.0018 
0.0809 ± 0.0031 



0.137 



+0.077 



-0.073 

0.255 (fixed) 
0.19 ± 0.11 

0.338 (fixed) 
0.25 ± 0.15 
0.322 (fixed) 
0.486 (fixed) 
0.289 (fixed) 



Tabic 5. Physical parameters 



Parameter 


D13 


This work (w/o b prior) 


This work (w/ 6 prior) 


Stellar parameters 
Stellar mass, Ms [M©] 
Stellar radius, Rs [R©] 
Stellar density, ps [p©] 
Stellar surface gravity, logps [cgs] 


n e:oq+0.047 

0.568^^31 

2.9it^;l 

4.658 ± 0.030 


n rc:7+0.028 
U.OO(_o.o20 

0.526«:»?? 

q oq+0.14 
0.00_o 32 

4 741+0.009 

'I- '^--0.019 


n c;q4+o.029 

n 56^+0 024 

u.ooo_o 020 

q q9 + 0.24 
'5"J^-0.27 

4 710+0016 

'±./lU_o 019 


Planetary parameters 
Orbital period, P [days] 
Semi-major axis, a [AU] 
Planetary mass, Mp [M0] 
Planetary radius, Rp [Re] 
Planetary density, pp [g cm~^] 
Planetary surface gravity, gp [m s~^] 


3.33671 ± 0.00005 
n n^^c;7+o-ooo96 

U.UOJOI -0.00100 

13.9llj 
4 83+°-22 

'i.OO_0 21 

72+0 " 
'J- '^-0.12 

r 7ir+0.85 


3.336652 ± 0.000005 

U.UaoyO_o. 00044 

14.1 ± 1.3 

'i-'J^-O.lO 

0.94 ± 0.12 

7 9i;+0.75 


3.336651 ± 0.000005 

U.UOT/4_0. 00054 

14.6 ± 1.4 

0.80 ± 0.11 
6.58lo;70 



